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Abstract. We investigate the electromagnetic response of a pair of complementary 
bi-anisotropic media, which consist of a medium with positive refractive index (+£, 
+fi, +£) and a medium with negative refractive index(— e, —fj,, — £). We show that 
this idealized system has peculiar imaging properties in that it reproduces images of 
a source, in principle, with unlimited resolution. We then consider an infinite array of 
line sources regularly spaced in a one-dimensional photonic crystal (PC) consisting of 
2n-layers of bi-anisotropic complementary media. Using coordinate transformation, 
we map this system into 2D corner chiral lenses of 2n heterogeneous anisotropic 
complementary media sharing a vertex, within which light circles around closed 
trajectories. Alternatively, one can consider corner lenses with homogeneous isotropic 
media and map them onto one dimensional PCs with heterogeneous bi-anisotropic 
layers. Interestingly, such complementary media are described by scalar, or matrix 
valued, sign-shifting parameters, which satisfy a generalized lens theorem, which can 
be derived using Fourier series solutions of the Maxwell's equations (in the former 
case), or from space-time symmetry arguments (in the latter case). Also of interest 
are 2D periodic checkerboards alternating rectangular cells of complementary media 
which are such that one point source in one cell gives rise to an infinite set of images 
with an image in every other cell. Such checkerboards can be mapped into a class of 
three-dimensional corner lenses of complementary bi-anisotropic media. 
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Figure 1. (a) Configuration for the focusing of a chiral slab for an object in air. 
Subscripts '+' and '-' of the electric field spectra indicate the left- and right-going 
wave streams, respectively, (b) Plot of \/e 2 + h 2 for a time harmonic electric line source 
radiating in the chiral slab by Finite element method. 



1. Introduction 

Jin and He recently proposed to use a slab of chiral medium in order to achieve a 
refocussing effect via negative refraction pQ. These authors noted that at the interface 
between a conventional medium with permittivity Eq and permeability no, and a chiral 
medium with permittivity E\ and permeability yU 1; and chirality k, the refractive indices 
for circularly polarized waves associated with wavenumbers k± = u(Ky/fioeo ± y/fiiEi) 
are n± = y (/ii^i) / (/io^o) ± K - Indeed, when a wave is obliquely incident upon such 
an interface, and if n^eofio > Js\iL\, then n_ is clearly negative, so that circularly 
polarized waves with wavenumber k_ undergo some form of negative refraction, as shown 
in figure [l] (rays in red color), which is the earlier result of Jin and He[lJ. Wiltshire and 
Pendry also proved a strong chiral in chiral Swiss rolls, which can be larger than either 
the permittivity or the permeability [2J. However, it is also legitimate to ask whether 
one can achieve similar focussing effects with more general bianisotropic media with 
sign shifting in permittivity, permeability and bianisotropy. Furthermore, one can also 
consider heterogeneous anisotropic parameters. One way to answer such questions is to 
consider complementary media which cancel out optical space [3], and to extend their 
analysis to chiral complementary media. Using the powerful tools of transformation 
optics, one can also extend the design of generalized perfect lenses and corners to chiral 
and bianisotropic media. Finite element methods shall be used to numerically explore 
such optical systems. These analytical solutions of the singular systems considered by us 
represent important benchmarks for the validation of non-trivial numerical calculations. 



2. A Generalized Lens Theorem for chiral media 

2.1. Proof via Fourier analysis 

The original 'perfect lens' as introduced by Victor Veselago[1968] [5] is defined by 
a slab of material with e — — 1, jj, — — 1, however, a much more general condition as 
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Figure 2. A complementary bianisotropic medium. 



proposed by Pendry and Ramakrishna in 2003 [3] is the system described by, 
E\ = +e(x,y), = +fi(x,y), — d < z < 

e 2 = -e(x, y), /i 2 = -fi>(x, y), < z < d 1 ' ' ' 

Within this framework, a new theorem can be stated to prove the perfect lens: it is 
enough to point out that two complementary media have optical sum zero (the media 
in the region —d < z < d behave as though they had zero thickness). 

However, when there is chirality in the complementary media, do they still work as 
a perfect lens? A new generalized lens theorem should be considered. In what follows, 
we provide a positive answer for this kind of complementary bi-anisotropic media, 

A system with two complementary bianisotropic media is shown in figure |2j 

ei = +e(ar,y), fix = +n(x,y), & = +£(x, y), -d<z<0 
e 2 = - £ (x, y), [i 2 = -fi(x, y), 6 = -£0, v), < z <d 



(2) 



Let us assume constitutive relations in a bianisotropic medium as follows, 

D = eE + z£H, B = /jlH. — i£E (3) 

with £ the chirality parameter, instead of k in He pQ. 

Maxwell's equations take the form (known as Maxwell- Tellegen's equations), 
„ „ <9H dE 

„ dE .m 

Let us decompose the electric fields in Fourier series, 
Ex(x, y, z) — exp(+ik\ z z) J] Ei(A; :c , k y ) exp(ik x x + ik y y), —d<z<0 

fax i^y 

E 2 (x,y, z) = exp(+ik 2z z) J2 E 2 (k x ,k y ) exp(ik x x + ik y y), 0<z<d 



(4) 



(5) 



The same holds true for the magnetic fields. These fields should satisfy some boundary 
conditions: The tangential components of both electric and magnetic fields are 
continuous across the interface at z = 0. 



Focussing light in a bi- anisotropic slab with negatively refracting materials 4 

Plugging the Fourier series ^ for the electric and magnetic fields into Q, and 
using the frequency domain notation d/dt = iu, we get, 

k u z x [E lx (k x , k y )x + Ei y (k x , k y )y] + (k x x + k y y) x E lz (k x , k y )z 

tufa, K> K, K) £ + H M, KM 



^ £,i(.kx,k y ,k x ,k y )[Ei x (k x ,k y )x + E ly (k x ,k y )y] 

b' b> 

n, x ,Ky 

(k x x + k y y) x [E lx (k x , k y )x + E ly (k x , k y )y] 

= us ^ Hi{k x , k y , k x , k y )Hi z (k x , k y , k x , k y )z 

— ico 2^ £,i{kx, k y , k x , k y )Ei z (k x , k y , k x , k y )z 

™x > K y 

k\ z z x \Hi x (k x , ky)x + H\y(k x , k y )y] + ik x x + k y y) x H\ z (k x , k y )z 
= ~ £ £ i(h, hi K> tf v )[Eix(K, k 'y)x + E iy(K> k' y )y] 

b 1 b 1 
b' h> 

(k x x + k y y) x [H lx (k Xj k y )x + H ly (k x , k y )y] 

e i{k>xi kyi k XJ k y )Hx z (k X ) kyj k x7 ky)z 



b> b> 

K x ,&y 



ico 2^ ky^k XJ k y )Hi z {k x ^ ky^ k x j k y )z 



b' b 1 



(6) 



If we now substitute, 

E2x{k x , k y ) = Ex x [k x , ky), E2 y {k x , k y ) = Ei y [k x , k y ), E2 Z [k x , k y ) = —E\ z ik x , k y 

H2x{k x ,ky) = Hi x (k x ,ky), H2 y {k x ,k y ) = H\ y (k x ,k y ), H2 Z {k x ,k y ) = H\ z {k x ,k y ), 

£2{k x , k y , k x , k y ) = —£i(k x , k y , k x , k y ), ^2(k x , k y , k x , k y ) = —fj,i(k x , k y , k x , k y ), 

£ 2 (k x , k y , k x , k y ) = —^i(k x , k y , k x , k y ), k2 Z = —k\ z (7) 
into (|6]), we obtain, 

- k 2z z x [E 2x (k x , k y )x + E 2y (k x , k y )y) - (k x x + k y y) x E 2z (k x , k y )z 

= -"E M k *> h> Ki k' y )[H 2x (k' x , k' y )x + H 2 y{k' x , k'y)y] 

b* b' 

+ Uh, k y , k' x , k' y )[E 2x (k' x , k' y )x + E 2y (k' x , k' y )y] 

h> b 1 

,Ky 

(k x x + k y y) x [E 2x {k x , k y )x + E 2y {k x , k y )y] 



ui ^ fi2{k x , k y , k x , k y )H2 Z (k x , k y , k x , k y )z 
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— %uj y j £2(^5 kyi k XJ ky)E2z(k X ) k y , k XJ k y )z 

b 1 b 1 

-k 2z z x [H 2x (k x , k y )x + H 2y (k x , k y )y] - (k x x + k y y) x H 2z (k x , k y )z 
= u; £ k y , K, k' y )[E 2x (k' x , k' y )x + E 2y (k' x , k' y )y] 

b f b' 

+ iu J2 6(fcx> k v , K, k y )[H 2x (k x , k' y )x + H 2y (k x , k' y )y] 

k' x ,fCy 

(k x x + k y y) x [H 2x (k x , k y )x + H 2y (k XJ k y )y] 

— ^2(^x5 ky-> k XJ ky)H2z{kxj k y , k XJ k y )z 



b 1 b' 



-i^Yj &(k x , ky, k' x , k' y )H 2z (k x , k y , k' x , k' y )z 

b 1 b' 



These equations can be rearranged as, 

k 2z z x [E 2x (k x , k y )x + E 2y (k x , k y )y] + (k x x + k y y) x E 2z (k x , k y )z 

= n 2 {k x ,k y ,k x ,k' y )[H 2x (k' x ,k' y )x + H 2y {k' x , k' y )y] 



b' bf 



k y , k XJ ky)[E 2x (k XJ k' y )x + E 2y (k x , k l y )y] 

b' b* 
n, x ,ruy 

(k x x + k y y) x [E 2x (k x , k y )x + E 2y (k x , k y )y] 

= u foiK, k y , k x , k' y )H 2z (k x , k y , k' x , k' y )z 

b' b' 

l^ X ^y 

-w &(k x ,k y ,k x ,k' y )E 2z (k x ,ky,k x ,k y )z 

k' x ,ky 

k 2z i x [H 2x (k x , k y )x + H 2y (k x , k y )y) + (k x x + k y y) x H 2z (k x , k y )z 
• £ 2<>x, ky, k' x , k' y ) [E 2x (k' x , k' y )x + E 2y (k' x , k y )y] 



—uj 
iuj 



&(k x , k y , k' x , k' y )[H 2x (k' x , k' y )x + H 2y (k' x , k' y )y] 

b' b* 

n, x ,Ky 

{k x x + kyij) x \H 2x {k x , ky )x + H 2y (k x , k y )y] 

= -u e i(k x , k y , k' x , k' y )H 2z (k x , k y , k' x , k' y )z 

£,2(k x , ky, k x , k y )H 2z {k x , k y , k x , k y 



IUO y t2[K x ,K v ,K„,K„,)tt 2 AK x ,K„,K„,K y )Z 

b' b* 



(8) 



(9) 

Equations ^ and ^ hold for the fields in chiral media with parameters and 
£ 2,^2,£,2 respectively. The fields in the two regions match across the interface due to 
the continuous boundary condition, while there is a sign-shifting between components 
along z direction, hence, 

E(x, y,z = +d) = E(x, y) exp (—ik z d) 
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E(x, y) exp (ik z (-d)) = E(x, y, z = -d), d > 0. (10) 



( 10 ) holds for each k z and hence the sums over k z , the same result holds for the magnetic 
fields. We would like now to extend these results to anisotropic tensors of permittivity, 
permeability and chirality. This can be done using the same kind of techniques as above, 
but computations become tedious. We therefore propose another, more subtle, but more 
concise, derivation of the generalized lens theorem using space-time symmetries. 

2.2. Proof via space-time symmetries 

Rewrite Maxwell- Tellegen's equations Q of a bianisotropic medium as, 
„ „ d(£E) d(uH) 

Wt '*^ (U) 
at at 



Let us list some transformations which leave the equation (11 ) invariant (see [6] for the 
case when £ = 0): 

51 (Generalized conformal invariance): E — > AE, H — > AH, ^T 1 — > A[i~ 1 A~ 1 , 
e — > AeA^ 1 , £ — > A^A^ 1 , where A is invertible and an element of GL 3 (R) (a group of 
3x3 linear operators). 

52 (Generalized duality): E — > — H, H — > — E, (iff /x = e). 

53 (Parity invariance): r — > — r, where r = [x, y, z], Maxwell- Tellegen's equations in 



(11) turn out to be, 



-V x E + i- 



dt dt /-| 2 x 



in order to preserve equation (11), E — > — E, H — > — H, e — > — e, /i — > — /i, £ — > — £ 
should be chosen to cancel out the minus introduced by — r. 

54 (Time reversal): t — > —t, then E — > — E (or H — > — H) and £ — > — £ (similar idea as 
Parity invariance). 

55 Any additional space-time symmetries. 

The combination of any of these symmetries is again a symmetry of the system of 
equations. In other words, if the fields in a particular region of space can be mapped on 
to another region of space through the symmetry transformations S1-S6 while preserving 
the respective boundary conditions, then the transformed fields solve the field equations 
whenever the original fields do. This law will also be proved in the following part with 
coordinates transformation. 

Let us consider a homogeneous slab in figure [2] with permittivity, permeability and 
chirality tensors in the region of — d < z < 0, 



Vxx Vxy V X z 



Vl = Vy X Vyy Vy Z , V = £, fj,, £ (13) 



V Z x ^zy V z 
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Figure 3. An alternative pair of complementary chiral media, each cancelling the 
effect of the other, that means complementary halves sum to zero. 



For a source at the interface and propagation along the z direction, we use the symmetry 
operations S3, followed by SI with, 

-10 

(14) 



A 














1 



We will call this sequence of operations a mirror operation. The choice of A preserves 
the continuity of En and Hn across the boundary. Then the complementary medium of 
< z < d is, 





^xx 


^xy "xz 






v 2 = 


Vyx 


— Vyy Vy Z 


, V = £, fl, £ 


(15) 




^zx 


Vzy V Z z 







Similar to the complementary theory proposed by Pendry and Ramakrishna, a schematic 
diagram for the bianisotropic complementary media is indicated in figure [3| two 
complementary bianisotropic media have an optical sum of zero. The optical response 
on the two sides of a structure with complementary media can be calculated by cutting 
out the complementary media and closing the gap between the two sides. 

The two derivations of the generalized lens theorem are underpinned by the 
invariance of the Maxwell's equations under geometric changes. If one wants to study 
other classes of generalized lenses, say in polar or spherical coordinates, some more 
elaborate mathematical tools are needed, as we shall now see. 



3. Transformation optics applied to the design of corner lenses for 
bianisotropic media 

3.1. General Coordinates Transformation 

Considering a new coordinate system described by[l], 

Qi(x,y,z), q 2 (x,y,z), q 3 (x,y,z) (16) 

Lines of constant g 2 , q$ define the generalized qi axes, and so on. It should be noted 
that we can produce a desired mesh by choosing the proper coordinates transformation. 
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Figure 4. A simple cubic lattice in one coordinate system (left) maps into a distorted 
mesh in the other coordinate system(right). 



For example, if we define a set of points by equal increments along the new coordinates 
qi, q2 and q 3 axes, these will make the mesh in the original coordinates system look like 
distorted as shown in figure |4j 

Since Maxwell- Tellegen's equations have the forms of Q in a Cartesian coordinate 
system, while for a new coordinate system of q\, qi, q 3 (which can be non-orthogonal), 
one may wonder what will happen to Maxwell- Tellegen's equations, i.e. will they preserve 
their form like, 



V,xE 



V,xH 



<9H 

dt ~ 
dE 



it 



dE 



dt' 



(17) 



dt * dt 

where E, H are re-normalized electric and magnetic fields, and e, fi and £ are in general 
tensors? While this has been discussed by Ward and Pendry for anisotropic systems [I] , 
we present here a generalization of those results to fully bianisotropic systems. 

The proof of the form invariance of the Maxwell- Tellegen's equations is detailed 
below: 

First, we define three units vectors. Let u 1; u 2 , u 3 be along the axes of qi, q 2 , q 3 . 
The transformation relation between [x,y,z] and [q±, q 2 , q 3 ] is obviously, 



(18) 





dx 


dx 


dx 




dx 




dqi 


dq 2 


dq 3 




dqi 


dy 




dy 


dy 


dy 








dqi 


dq 2 


dq 3 




dq 2 


dz 




dz 


dz 


dz 




dq 3 




. dqi 


dq 2 


dq 3 . 





the length of a line element in the new coordinate system can be derived from 

ds 2 = dx 2 + dy 2 + dz 2 = Qndq 2 + Q 22 dq\ + Q 33 dq\ 

+ 2Qi 2 dqidq 2 + 2Qi 3 dqidq 3 + 2Q 23 dq 2 dq 3 (19) 

with 



^ dx dx dy dy dz dz 

Qi ■ = 1 — — H 

%0 dqidqj dqidqj dq { dqj 



(20) 
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E l Q x dq l u l + dq 2 —(E l Q ] dq l ) 



Q 3 dq 3 uJ /Q 2 dq 2 u 2 



E 2 Q 2 dq 2 Uj. 



E,Q,dq 2 u 2 + dq l - — (E 2 Q,dq 2 ) 



E l Q 1 dq 1 u, 



Figure 5. (a) Small element resembling a parallelepiped, (b) Integration path for 
finding V x E. 



Let us now note that the electric field E and magnetic field H can be expressed in 
terms of contravariant components, 

U = U 1 ^ + U 2 u 2 + U 3 u 3l U = E, H (21) 

and conversely, they can be expressed in terms of the covariant components through 





' u 1 ' 




g- 1 


u 2 






u 3 





U X ■ U X U X ■ U 2 U! • U 3 
U 2 ■ Ui U 2 ■ U 2 U 2 • U 3 

U 3 ■ Ux u 3 ■ u 2 u 3 • u 3 

with g a 3*3 matrix and 

U x = U • ui, U 2 = U • u a , U 3 = U • u 3 . 





" U 1 ' 




' u x ' 




U 2 




u 2 




U* 




U 3 



(22) 



(23) 



according to ( 22 ) , 



U l = Y,9 lJ U 3 (24) 
i=i 

We take a small element which resembles a parallelepiped (left panel in figure [5]), 
and consider the projection of V x E onto the normal to the u x — u 2 plane by taking a 
line integral along the ui — u 2 parallelogram (right panel of figure p| and invoking 
Stokes' theorem, 

" V x E • dS = /e • dl (25) 
s Ji 

where S is the surface of Ui — u 2 plane and I is the boundary around surface S. 
Consequently, 

d d 
(yxE)-( y u 1 xu 2 )Q 1 dq 1 Q 2 dq 2 = dq l — (E 2 Q 2 dq 2 )-dq 2 — (E l Q l dq l )(2Q) 

c% dq 2 

on the other hand, the normalized electric field in new coordinate system is defined as, 
Ei = QiEi, (2 = 1,2,3) (27) 



Hence equation (26) can be simplified as, 



(V x E) ■ ( Ul x u^Q.Q, = 7 r i - T i = (V,xE) 3 

dqi dq 2 



(2* 
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the term in the right-hand side of this equation is the third component of curl in the 
new coordinate system. Then applying Q to V x E in the left-hand side of equation 

(28) we obtain 

OH <9E 

(V x E) • (m x u 2 )QiQ 2 = ■ (ui x u 2 )Q 1 Q 2 -i^—- (m x u 2 )Q!Q 2 (29) 

Here it should be pointed out that the equation (15) in the paper of Ward and 
Pendry has a flaw for non-orthogonal coordinates system q\,q 2 ,qz- However, this does 
affect their result. Let us now substitute (21) combining with (24) for both E and H to 

(29) and obtain 



<9H <9E 

(V x E) ■ (ui x u 2 )QiQ 2 = H-^r ■ ( u i x u 2)<2i<?2 - i^-^r ■ ( u i x u 2 )QiQ 2 



3 QJJ ■ 

3=1 



3 BE 

^E^ 3i ^/ u 3-(uixu 2 )g 1 g 2 (30) 



3=1 



ot 



define 



and 



vg K \vix ■ (u 2 x u 3 )\QiQ 2 Q3(QiQj) \ 



U; 



'3 QjUj 

hence (30) turns out to be 



(VxE)-(u 1 xu 2 )Q 1 Q 2 = ^/i 

3=1 



dt 



3=1 



dt 



comparing with (28), we finally obtain 

<9H . t dE 



Vg X E 



dt ^ dt ' 

Similarly, due to the symmetry between E and H fields, we have 



V,xH 



: ~dt 



it 



~dt 



(31) 
(32) 

(33) 

(34) 
(35) 



while e %] also satisfies (31). 

These results indicate that the form of Maxwell- Tellegen's equations is preserved 
under a coordinates transformation, while the definitions of e, fi and £ are simply 
changed. In other words, no matter what is the coordinate system, we can always analyze 
the electromagnetic property of a structure through the Maxwell- Tellgen's equations 
with redefined parameters. 



3.2. ID PCs to 2D corner lenses 

In the previous section, we have analyzed a system with two complementary regions: 
the optical changes through these two regions is zero due to the fact that two regions 
are inverted mirror images of each other. In the following sections, we would like to 
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(a) (b) (c) ' 

Figure 6. (a) A pair of 2D corners of complementary chiral media can focus a source 
back on to itself, (b) The cylindrical coordinates, (c) A layered system with a periodic 
set of sources. 



give a further generalization of the cancelation principle combining with coordinates 
transformation. 

Using a general method of coordinates mapping, one can map Maxwell- Tellgen's 
equations to other geometries and obtain a class of generalized perfect lenses in 
curvilinear coordinates. For instance, two negative two-dimensional corners, within 
which light radiating from a line source is bent around a closed trajectory and is 
refocused back onto itself [9], can be mapped to a layered structure consisting of layers 
of complementary bianisotropic media together with a periodic set of line sources, as 
shown in figure |6j Alternatively, the inverse transformation from a layered structure to 
2D corners is available. 

Here, the parameters of the chiral medium have been considered in the most general 
case which can take any values. Specifically, permittivity and permeability of a perfect 
lens proposed by Pendry is supposed to be -1, while for a perfect lens with chiral medium, 



its refractive index is n = \Jejj/eoHo ± co£, that Co is the speed of light in vacuum and 
SqUq = Cq 2 , hence by taking proper values of the permittivity, permeability and chirality, 
we can obtain n, correspondingly, its complementary medium has — n, for a source placed 
in any one cell of the 2D corners, we have an image point in every cell of the corners. 

Here, we would like to give a proof of the mapping from 2D corners onto periodic 
layered systems. For this, let us consider, 

e(ip) = +e(<p), fi{(p) = +n{<p), £(<p) = +£(<P), ?r/2 < y? < vr 

e(<p) = -e((p), fi((p) = -fi((p), £(¥>) = -£(¥>), < <p < vr/2 

then let us introduce a mapping of coordinates that takes structure in the left panel of 
figure [6] into the structure of the right panel. 

x = r cos ip e l / l ° , y = r sin ip e l ^ l ° , z = Z (37) 

while r , <p and Z are the cylindrical coordinates, in which the layered structure is 
defined. In this frame, 

^ii 
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where, 

Qi = r /lo\Je 2l / l ° cos 2 <p + e 2l / l ° sin 2 ip = r /l Q e l/l ° 



Q<p — r y e 2l / l ° cos 2 ip + e 2l / l ° sin ip = r e 

Qz = i 

QiQ v Qz = r 2 /l e 1 ^ 



l/lo 



(39) 



and hence, 



(41) 



ei = l si, e v = l 1 s (p , e z = r 2 /l e 2l/lo e z 

p,i = lolJ,i, p, v = l 1 fi v , jx z = r 2 /l e 2l/lo fi Z (40) 
ii = Ui, i? = la%, lz = rl/l Q e 2l ' l °i z 
substituting from ( |3~6| ) and setting l = 1 gives, 
7r/2 < (p < 7i : 

i l = E lfi = e(ip), e z = r^e 2l E((p) 
h = K = Mv)> frz = r%e 2l n{p) 

ii = i<P = m, iz = r 2 e 21 ^) 

< p < tt/2 : 

ei = e v = -e(<p), e z = -r 2 e 2l E((p) 
h = K = = -r 2 e 2l fi{<p) 

1 = Z v = -Z{<p), lz = -rle 2l i{p) 

which can be indicated in the ID periodic layers with an alternation of complementary 
bianisotropic medium. 

3.3. From 2D checkerboards to 3D corner lenses 

Figure [7] shows a 2D checkerboard consisting an alteration of rectangular blocks 
of complementary bianisotropic media. The axes along the checkerboard directions are 
defined by 9 and 0, along which the periodic of the checkerboard is tt, the direction 
normal to this plane is denoted by I. The structure is defined, 

e(9, 0) = +e(6, 0), »(6, 0) = 0), f (0, 0) = +£(6, 0), V (9, 0) e positive, 

e(0,0) = -£(0,0), /i(0,0) = -/i(0,0), £(0,0) = -£(0,0), V(0,0) E negative, 

where 

positive = (0, tt/2) x (0, tt/2) U (tt/2, tt) x (tt/2, tt), 

(43) 

negative = (0, tt/2) x (tt/2, tt) U (tt/2, tt) x (0, tt/2), 

Obviously, a refocussing can be observed in each cell while the rays starting from a 
single source in a positive cell as shown in figure [7j This phenomenon can be unveiled 
by the generalized lens theorem, let us consider the imaging along the 9 direction, in 
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Figure 7. An infinite checkerboard with optically complementary cells. When there is 
a source in one cell of the checkerboard, then according to the ray routes, images are 
recreated in every other cell. 




(a) (b) (°) 

Figure 8. (a) 2D checkerboard with optically complementary cells, (b) The spherical 
coordinates. (c) Cubic corner lens. 



which direction the condition of complementarity is satisfied for the layers, hence, for a 
source at 9 — 9i, <fi — 4>x in a positive cell, we can obtain a set of images along <fi = <fii 
line at 9 = ±(nii + 9i), with n a positive integer. Similarly, a set of images along </> 
direction can be obtained, and finally an image in every cell of the checkerboard can be 
achieved corresponding to the source placed in any one cell. 

According to proper coordinates transformation, a 2D corner can be mapped to a 
ID lens which satisfied the complementary conditions. Furthermore, let us consider the 
mapping algorithm from a 2D checkerboard to 3D cubic corner as shown in figure [8] As 
shown in the left panel of figure [8j the checkerboard consisting of a doubly periodic set 
of point sources with period 27r along and period 7r along 9, the coordinates of which 
is denoted by [I, the right panel is corresponding to a 3D cubic corner extending 
to infinity, of which the coordinates are denoted by Cartesian coordinates [xi,x%, X3]. 
According to the spherical coordinate system indicated in the middle panel, the two 
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coordinates satisfy. 

I = 
9 = 



X3 



arccos 



arctan 



/->-» 2 | ,-y-i 2 I y ■ 2 



), 



(44) 



where ro is a scale factor and can be considered as the radial position of the source. I 
denotes the radial (logarithmic) coordinate and 9, are the longitudinal and azimuthal 
coordinates; we can then generate the corresponding cubic corner system in [x\, x%, X3]. 

Considering a checkerboard with homogeneous and isotropic cells, then according 
to the coordinates transformation algorithm, we have 



/ 2 <y. 4 ™ 2 

l o' X 1 


+ 




^2 ^,2 


+ 


rp 2 

x 2 




r 4 


Cac? + A) 


(xf + x 2 ,) 2 


72^,4,^.2 
f 0' X 2 


+ 




, >-» 2 , y , 2 


+ 


x 2 




r 4 


(x? + x|) 


(x 2 + x 2 ) 2 


72^4^2 
l 0' X 3 


+ 


7 ■ 2 


+ x 2 






r 4 




r 4 ' 







Qx = 

Q2 = 

Q 3 = 

with r = ^x 2 + x 2 . + x 2 . Assuming l = 1, r = 1 and 

Xi = r sin# cos0, X2 = r sin 9 sin 0, X3 
(45) can be recasted as, 

1 / „ , sin 2 



(45) 



r cos 9, 



Qi 



< cos^ 



+ 



sin 2 9 ' 



ft 



' sin 



+ 



cos 



sin 2 9 



Q-3 



(46) 



(47) 



on the other hand, we have (38) as the new set of the parameters, hence let us substitute 
(47) into (38), we can deduce 

Vi sin 2 9 sin 2 + cos 2 



Vi 

v 3 



r sin 2 9 cos 2 + sin 2 
v 2 , sin 2 9 cos 2 + sin 2 
r sin 2 # sin 2 + cos 2 

^3 



(4£ 



r sin 6 1 



sin 2 cos 2 + sin 2 0\/ sin 2 9 sin 2 + cos 2 



where v = e, fx, £. It is noted that the transverse r is irrelevant to the imaging. Now, 
choosing 



£1 
£3 



/*3 



6 = 
= /i2 

6 = 



, sin 2 9 cos 2 + sin 2 < 
1 sin 2 9 sin 2 + cos 2 < 



= fh =6 
r sin 2 (9- 



(49) 



1 



1 



sin 2 # cos 2 + sin 2 a/cos 2 # sin 2 + cos 2 



then we obtain that (£1,82, £3), (^1,^2,^3) and (^1,^2,^3) to be identical to the 
checkerboard made of homogeneous, isotropic materials. So far, we have proved that 
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we can achieve a cubical corner lens made up of anisotropic, inhomogeneous materials 
by mapping from a homogeneous checkerboard. 



4. Numerical illustration of finite element method 



4-1- Diffraction problem for fully bianisotropic media 

In this section, we would like to give a numerical illustration to the structures 
consisting of complementary bianisotropic media through FEM (implemented in 
COMSOL MULTIPHYSICS). First, we would like to derive the diffraction equations 
for a bianisotropic medium, generally, we suppose those parameters in equation Q, 
such as permittivity, permeability and chirality are all tensors written as follows, 



^11 ^12 
^21 £22 








£33 



/ill 





1^12 
A*22 








1^33 



61 

61 




62 
£22 







£33 



(50) 



and also assume that E = (ei, e) ,H = (hi, h 2 , h) , which leads to 

V- {(^r + lre^- 1 t)~ 1 Ve}+iV- {(|r + |z k^ 1 /^V/i} 

= u 2 e 33 e + iu 2 £ 33 h 
*V ■ {(^ T + It 1 e^-'Ve} + V ■ {(^ + |r ^r" 1 t^^Vh} 

= iu 2 £ 33 e - u 2 fi 33 h 



(51) 



with 



e T 



£22 —£21 
—£12 £11 



-1-122 
Hl2 



fJ-21 
-/^ll 



T 



—&2 
62 



6l 
-6l 



(52) 



Based on this, we can set up any proposed model indicated by equation (51 ) in COMSOL 



MULTIPHYSICS by using PDE module. We note that such a model is required e.g. in 
order to implement perfectly matched layers (for which e, fi and £ are all tensors). In 
the sequel, we apply this numerical model to various types of complementary media. 



4-2. Complementary media with isotropic e, \i and £ 

In this part, we would like to give a numerical analysis to the structure of 
complementary bianisotropic media described by isotropic (sign-shifting) parameters 
e, n and £. First, we analyze an infinite chiral lens, which consists of the complementary 
bianisotropic medium as shown in the left panel of figure [9j while there is a line source 
lying in the top bianisotropic layer, in order to asymptotic the infinite half space, we have 
introduced the Floquet-Bloch boundary at the upper and lower interface, and perfect 
matched layers at the left and right sides of the structure. The FEM results is shown 
in the right panel, a refocussing of the source can be observed in each layer as indicted 
in previous theoretical analysis. Here, for the module in COMSOL, we have supposed 
e = e , /j, = /i and £/c = 0.99 respectively, with the normalized frequency of the source 
is / = 7 x 10 14 Hz. 
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Floquet Boundary 
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Floquet Boundary 



(c) 




Figure 9. Periodic ID complementary bianisotropic media: (a) Schematic view of the 
computational model with Perfectly Matched Layers (PMLs) on either sides of the 
structure, periodic conditions on top and bottom and a line source in the top layer; 
(b) 2D plot of \J e 2 + h 2 for a time harmonic electric line source radiating in a periodic 
set of complementary chiral lenses; (c) Plot of real part of electric field e; (d) Plot of 
real part of magnetic field h. 



Furthermore, a 2D chiral corner is also considered as shown in figure 10 (the size 
of structure is in order of um), with a line source lying in one of corner with positive 
parameters, similarly, the source has been mapped to the other corners as an image. The 
parameters of the media are supposed to be e = Eq, fi = p,o and £/co = 0.99 respectively, 
where the normalized frequency of the source is / = 6 x 10 14 Hz. 

In the same way, a triangle corner lens as shown in figure [TT] is simulated, with 
e = So, /i = /io and £/co = 0.99 respectively, and the normalized frequency of the source 
is / = 7 x 10 14 Hz. 

At last, a checkerboard and its numerical results by COMSOL are shown in figure 
T2| the permittivity and permeability are e = e , \i = /x , while the chirality is 
£/c = 0.99, the frequency is / = 3 x 10 14 Hz. 

One should note that absence of absorption in sign-shifting bianisotropic media 
means the Maxwell- Tellegen operator is no longer elliptic and hence Lax-Milgrams 
lemma ensuring existence and uniqueness of the solution to the diffraction problem 
is no longer applicable. Nevertheless, the Fredholm alternative ensures that if there is a 
solution for any applied electromagnetic field whose tangential component is continuous 
through interfaces separating complementary media (hence exhibiting two anti-parallel 
wave-vectors at both sides of such interfaces) and satisfying periodic conditions on the 



Focussing light in a bi- anisotropic slab with negatively refracting materials 



17 





PML 








• 




: 

1 

/// 




-s,-ft,-4 






— 




PML 





(a) 



(c) 













H 






n 













I 



0.06 
0.W 
0,02 



(b) 



I 



™0 

I 



i 



(d) 



Figure 10. 2D complementary bianisotropic corner: (a) Schematic view of the 
computational model with Perfectly Matched Layers (PMLs) on either sides of the 
structure, and a line source in one of the corner; (b) 2D plot of Ve 2 + h 2 for a time 
harmonic electric line source radiating in a 2D complementary bianisotropic corners; (c) 
Plot of real part of electric field e; (d) Plot of real part of magnetic field h. 



unit cell edges, it will be an eigenf unction for the problem as shown in figure |12[b)- (d) . 
The convergence of the Finite Element algorithm depends crucially upon the regularity 
of the boundaries between complementary media with some ill-posedness occurring 
in certain corner cases for a refractive index contrast of —1 when the coercive plus 
compact mathematical framework breaks down, see [TU] for the case when £ = 0. Thus, 
our computer calculations represent a benchmark for validation of numerical results 
involving bianisotropic materials in very singular conditions such as those presented 
here. 

4-3. Numerical results for anisotropic e, /x and £ 

Let us now look at the case of a cylindrical lens with tensors of permittivity, 
permeability and chirality as follows 

V! = diagfv^, #, vjp], v = e,fi,( (53) 
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Figure 11. 2D complementary bianisotropic corner with an alteration of triangle 
shape: (a) Schematic view of the computational model with Perfectly Matched Layers 
(PMLs) on either sides of the structure, and bottom and a line source in the top layer; 
(b) 2D plot of Ve 2 + h 2 for a time harmonic electric line source radiating in a periodic 
set of complementary bianisotropic corners; (c) Plot of real part of electric field e; (d) 
Plot of real part of magnetic field h. 



with the index 7 = 1,2,3 corresponding to the three regions of the cylindrical lens 



+1, tiW = +1, «w = +1, r < n 
— 1, v^p = —1, = —r\jr^, t*i < r < r 2 
+ 1, v^' = +1, = rf/rf, r 2 < r < r 3 

where region 3 is the vacuum surrounding the cylindrical lens, region 2 is filled 
with negative bianisotropic media, which is designed to optically cancel out region 3 
( r 3 — r \l T \)i and region 1 is interior part of lens; the whole structure acts that the 
region 1 is optically equivalent to vacuum of radius r 3 . Hence, for an object lying in 
region 1, its image will be magnified by a factor of r^/ri = r\jr\ in region 3. 

It should be noted that when the source placed at r > rf/r 2 , images will be formed 
in both the annulus and the outer vacuum space as shown in figure fl3tb); otherwise, 



;(2) 
X 

,(3) 



(1) 



(54) 



when r < r^/r 2 , only an image will be formed in region 3 in figure 13 (c); while when a 
source is placed in the annulus, then there will be two images respectively in region 1 
and 3 as in figure 13 (d). where frequency / = 8.7 x 10 14 Hz, and the radius r\ = 0.2um, 
r 2 = OAum. 
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Figure 12. An infinite checkerboard consisting with complementary bianisotropic 
media: (a) Schematic view of the computational model with Floquet Boundary 
conditions in either side of the unit cell; (b)-(d) Plots of y/e 2 + h 2 for a line source 
radiating in an infinite checkerboard of complementary bianisotropic media. 



It is interesting to note that anomalous resonances on either sides of the cylindrical 
lens crucially depend upon absorption within the shell, see [IT] , and we have checked 
that introducing some small imaginary part to one component of the permittivity, 
permeability or chiral tensor is enough to remove the singularities on the shell's 
boundaries. These anomalous resonances can be used in a way similar to what was 
done in [12J to cloak some dipole sources outside a shell with negative permittivty, but 
this lies beyond the scope of the present paper. 



5. Concluding remarks 

In conclusion, we have derived a generalized perfect theorem that includes 
bianisotropic and chiral media in addition to the usual dielectric and magnetic media. We 
proposed two different methods to derive a generalized lens theorem for isotropic (section 
2.1) and anisotropic (section 2.2), sign-shifting permittivity, permeability and chirality 
admittance tensors. The geometric transformation technique (section 3.1) is generalized 
to include bianisotropic media as well and the two results are applied to analyze a 
variety of singular situations involving corners and wedges of such media. These results 
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Figure 13. A cylindrical lens consisting with complementary bianisotropic media: (a) 
Schematic view of the computational model with PML in either side of the unit cell; 
(b) Plots of \J e 2 + h 2 for a line source radiating in region 1 and r > rf/rs] (c)Plots of 



V e 2 + h 2 for a line source radiating in region 1 and r < r\/r2]{A) Plots of Ve 2 
a line source radiating in the annulus with r = r\ + [r-i 
axis from the origin. 



h 2 for 

ri)/2 along the horizontal 



are numerically simulated by finite element calculations using the COMSOL package 
and represent a benchmark for the validation of computations involving bianisotropic 
media in extremely singular circumstances. 
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